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ABSTRACT 



Aims. Stokes profiles emerging from the magnetized solar photosphere and observed by SOT/SP aboard the HINODE 
satellite present a variety of complex shapes. These are indicative of unresolved magnetic structures that have been 
overlooked in the inversion analyses performed so far. Here we present the first interpretation of the Stokes profile 
asymmetries measured in the Pel 630 nm lines by SOT/SP, in both quiet Sun internetwork (IN) and network regions. 
Methods. The inversion is carried out under the hypothesis of Micro-Structured Magnetized Atmosphere (MISMA), 
where the unresolved structure is assumed to be optically thin. We analyze a 29.52^^ x 31.70^^ subfield carefully selected 
to be representative of the properties of a 302^^ x 162^^ quiet Sun field-of-view at disk center. 

Results. The inversion code is able to reproduce the observed asymmetries in a very satisfactory way, including the 35 % 
of inverted profiles presenting large asymmetries. The inversion code interprets 25 % of inverted profiles as emerging 
from pixels in which both positive and negative polarities coexist. These pixels are located either in frontiers between 
opposite polarity patches or in very quiet regions. kG field strengths are found at the base of the photosphere in both 
network and IN; in the case of the latter, both kG fields and hG fields are admixed. When considering the magnetic 
properties at the mid photosphere most kG fields are gone, and the statistics is dominated by hG fields. According to 
the magnetic filling factors derived from the inversion, we constrain the magnetic field of only 4.5 % of the analyzed 
photosphere (and this percentage reduces to 1.3 % when referred to all pixels, including those with low polarization 
not analyzed). The rest of the plasma is consistent with the presence of weak fields not contributing to the detected 
polarization signals. The average flux densities derived in the full subfield and in IN regions are higher than the ones 
derived from the same dataset by Milne-Eddington inversion. 

Conclusions. The existence of large asymmetries in HINODE SOT/SP polarization profiles is uncovered. These are 
not negligible in quiet Sun data. The MISMA inversion code reproduces them in a satisfactory way, and provides a 
statistical description of the magnetized IN and network which partly differs and complements the results obtained so 
far. From this it follows the importance of having a complete interpretation of the line profile shapes. 

Key words. Sun: surface magnetism — Sun: photosphere — Techniques: polarimetric 



1. Introduction 



, The spectropolari meter SOT/SP (jLites et all 120011 : 
. iTsuneta et al.ll20Q8D ab oard the Japanese HINODE satel- 
' lite (|Kosugi et al.ll2007[ ) combines 0^32 angular resolution 
, and 10~^ polarimetric sensitivity to perform seeing-free 
■ full Stokes measurements of the polarized light emerging 
' from the magnetized solar photosphere. Since its launch 
in 2006, the SOT/SP instrument has allowed the solar 
community to investigate the solar surface magnetism 
under unprecedented stable conditions. 

HINODE SOT/SP observations have been largely ex- 
ploited to investigate the quiet Sun magnet ism. Its global 
propert ies have been descri bed by OjX)z cQ Suarez et al.l 
(|2QQ7aD, Lites et al. (2008), lAsensio RamosI (|2009[ ). and 
I Jin et al.l (|2009); while detailed analyses of its local prop- 
erties, in relation with the temporal evolution due to 
the interaction with plasma motions, h ave been carried 
out by [Ce nteno elal] (l2007f). [ Orozco Suare z et al. (2008), 
iNagata et a l. (2008), Shimi zu et all (|2008f ). iFischer et^ 
(|2009[ ) and lZhang et al.. i^2o6^ . 



The above studies have been performed exploiting in- 
version techniques for interpretation of the observed Stokes 
profiles. In most of the cases, the hypothesis of ME at- 
mosphere has been adopted to infer the properties of the 
magnetic field vector in HINODE resolution elements. It as- 
sumes that the polarization is produced in an atmosphere 
where the magnetic field vector is constant. This approx- 
imation is a good compromise to interpret large data sets 
in reasonable times, as the ones obtained by HINODE. 
However, it is important to realize that HINODE angular 
resolution is still too low to perform spectropolarimetric ob- 
servations of magnetic structures that can be rega rded as 
uniform (e.g., |SaiichezAlmeida'"2004"; "Stenfio"'2009). As we 
emphasize in § 13. If asymmetric Stokes profiles in HINODE 
quiet Sun observations are the rule rather than the ex- 
cept ion^_which_jni2l^ of unresolved struc- 
ture (jSanchez Almeida et al.lll996l , and references therein). 
The ME analyses cannot reproduce them, overlooking the 
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expected coexistence of several magnetic components in a 
single resolution ele ment Ej . 

Sanchez Almeida et alJ (|l996f ) put forward arguments 
for describing the atmosphere where the polarization 
is formed as a Micro- Structured Magnetic Atmosphere 
(MISMA), i.e., an atmosphere in which magnetic fields have 
structure smaller than the photon mean-free path at the so- 
lar photosphere 100 km). It can b e regarded as a limit- 
ing ca se of structures of all sizes (e.g. iLandi Degrinnocentil 
Il994al lbl). The MISMA hypothesis simplifies the radiative 
transfer, yet it provides realistic asymmetric Stokes pro- 
files. Under this hypothesis, the polarization from a sin- 
gle pixel in spectropolarimetric observations is equivalent 
to that produced by the average atmosphere. If one con- 
siders several components with diverse physical proper- 
ties (i.e., thermodynamics, plasma motions, and magnetic 
fields), the resulting spectrum is not a linear combination 
of Stokes profiles emerging from each magnetic compo- 
nent. Rather, the superposition is non-linear, giving rise 
to asymmetries with the properties observed in the Sun 
(e.g., with spectral lines that produce net circular polar- 
ization) . ISanchez Almeid a & Lites| (|2000i) showed how such 
a complex scenario can be properly adapted using a three 
component m odel MISMA, whose im plementation in an in- 
version code (S anchez Almeidalll997 ) allowed them to re- 
produce the full variety of profile asymmetries emerging 
from the quiet Sun when observed with the instrumenta- 
tion available at the time. Stokes profiles in IN and network 
observations performed with HINODE exhibit very impor- 
tant asymmetries which encouraged us to attempt the same 
MISMA analysis on these data. 

We present an orderly inversion of Stokes profiles ob- 
served by HINODE SOT/SP in IN and net work regions. It 
is performed using the inversion code in Sa nchez Almeidal 
(|l997). that has been recently employed to explain the 
reverse pol arity patches found by HINODE in sunspot 
penumbrae ([Sanchez Almeida fc Ichimotoll2009[ ). Our work 
represents the first analysis of quiet Sun SOT/SP data that 
is able to reproduce the asymmetries of the Stokes profiles. 
It allows us to obtain information contained in Stokes pro- 
files that are hidden to the ME analysis. In particular, we 
constrain the physical properties of the unresolved mag- 
netic structure creating the asymmetries. 

The paper is structured as follows: in § [2] we briefiy 
present the dataset and the selection of a sample subfield 
representative of the whole field-of-view (FOV). In §[3] the 
inversion hypotheses are exposed in detail, together with 
the adopted inversion strategy, and a few inversion exam- 
ples (§ 13. ip . The main results from the analysis are pre- 
sented in § m These are then discussed in § [5l An outhne 
of conclusions is given in § [6l 

2. Dataset and subfield selection 

We analyze full Stokes profiles of a 302^' x 162'^ portion 
of the solar photosphere observed at disk center on 2007 
March 10 between 11:37 and 14:34 UT. The spectropolari- 
metric measurements have been taken by the SOT/SP in- 
strument aboard HINODE in the two Fel 630 nm lines, 
with a wavelength sampling of 2.15 pm pixel" ^, and a spa- 

^ ME codes consider the existence of a field-free atmosphere 
portion i n the resolution element via a stray- li^ht filling facto r 
(see e.g. JSkumanich Litesll 19871 : lOrozco Suarez et al.ll2007al V 
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Fig. 1. Comparison between the statistical properties of the 
fuh 302'' X 162'' FOV (dashed line) and the statistical prop- 
erties of the selected 29.52" x 31.70" subfield (solid line). 
From top to bottom: histogram of circular polarization (C, 
first plot), total linear polarization (L, second plot), net- 
circular polarization {NCP^ third plot), and net-linear po- 
larization {NLP^ fourth plot). The error bars are defined 
as the square root of the number of counts in each bin. 



tial sampling of 0.1476" pixel" ^ and 0.1585" pixel" ^ along 
the east-west and south-north directio ns, resp ectively. The 
dataset has been alr eady analyzed by lOrozco S uarez e tldl 
(|2007aD . iLites et al.l ()2008.) and .Asensio Ramos. (;2009l ) to 
derive magnetic properties of IN and network regions. 

The data reduction and calibration have been per- 
form ed using the sp prep . pro routine available in solar 
soft (jlchimoto et al.l l2008f ) . For a correct absolute wave- 
length calibration, a correction for the gravitational redshift 
of 613 m s~^ has been taken into account. Using the po- 
larization signals in continuum wavelengths, we estimated 
a noise level of ay 1.1 x 10~^ Ic for Stokes V, and of 
(7q (7u — 1.2 X 10 ~^ Ic for Stokes Q and U {Ic stands for 
the continuum intensity). 

The inversion is carried out on a 29.52" x 31.70" subfield 
(equivalent to 200 x 200 pixels) of the fuh FOV. MISMA in- 
versions are time consuming, and this selection permits the 
in-depth analysis of a representative FOV in a reasonable 
timescale. The representative subfield has been selected 
based on two criteria. First, we consider invertible only 
those Stokes profiles presenting a maximum amplitude in 
Stokes V or Q or U above 4.5 times their noise. This check 
is performed in two different wavelength windows centered 
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Fig. 2. Main observational properties of the representative 29. 52^' x 31. 70^' subfield selected for in-depth analysis. 
Continuum intensity (upper- left panel), total linear polarization (L, upper-central panel), net-linear polarization {NLP, 
upper-right panel), COG magnetogram saturated at ±200 G (lower- left panel), total circular polarization (C, lower- 
central panel), and net-circular polarization {NCP, lower-right panel). The total polarization maps are both saturated 
at 1 pm, while the net-circular and net-linear polarization maps are saturated at 0.3. Black pixels in L, (7, NLP, and 
NCP maps represent regions with signals below the threshold for inversion. The bars on top represent the color palettes 
adopted for the two pairs of images below them 



on the two spectral linefl Second, the subfield must be a 
good representation of the whole set of invertible data at- 
tending to the statistical properties of the total circular po- 
larization {C = J Y^V(A)^dA//c), total linear polarization 
{L = J \/Q{X)'^ + /7(A)^dA//c), net-circular polarization 
{NCP = \AV\/ /|V(A)|dAj3, and net-linear polarization 
{NLP = {\AQ\ + \AU\)/ f \Q{X)\ + \U{X)\dX)^. The distri- 
butions employed to characterize the statistical properties 
of the full FOV and subfields use 30 bins to sample the full 
domain of each quantity. Bins with less than 16 counts in 
the case of the subfields, or 25 in the case of the full FOV, 
are discarded. In the subfield selection procedure, all pos- 
sible 29.52^' X 31.70'' subfields in the 302'' x 162" FOV are 
considered. The selected subfield is the one that minimizes 
the sum of differences between the subfield distributions 
and the ones defined over the 302" x 162" domain. Fig. [1] 
compares both the 302" x 162" statistics and the selected 
subfield statistics, while Fig. [2] contains the selected subfield 
continuum intensity represented tog ether with the center- 
of-gravity magnetogram (COG; see iRees fc Semell [1979 ). 
and the C, L, NCP, and NLP maps. 

3. Inversion hypothesis and strategy 

In order to assign physical properties to each point of the 
portion of photosphere under examination, we work out 



^ Each wavelength window is defined by two wavelengths 
(Aa and Xb) centered around the wavelength of the Stokes 
/ minimum (Am)- The width of the window centered on the 
Fel 630.15 nm line is 86 pm, while for Fei 630.25 nm the width 
is 64 pm. 

^ AS = f^"^ |'S'(A)|dA - r-^^ |6'(A)|dA, where S can be Q, U, 
or V, and Aa, As, and Am are defined in footnote ^. 



a model atmosphere producing Stokes profiles as close as 
possible to the observed ones. In the technical literature, 
this procedure is usually referred to as Stokes profile inver- 
sion. We perform the inversion under the MISMA hypoth- 
esis. In a Micro-Structured Magnetized Atmosphere, mag- 
netic fields vary over scales smaller that the photon mean- 
free-path at the solar photosp here (jSanchez Almeida et al.l 
'l996l; ^Sanchez Alm eidal 19971). Starting from this hypoth- 



esis, Sanchez Almeid a fc LitesI (|2000l ) derived by inversion 
eleven classes of MISMA models representative of typical 
profiles observed in IN and network regions at disk center. 

All these models are formed by three components: one of 
them is field-free, while the other two are magnetized. The 
field-free component represents the non-magnetized plasma 
in which the magnetic fields are embedded, while the two 
magnetized components allow us to model the coexistence, 
in the resolution element, of different magnetic fields con- 
tributing to the formation of the o bserved p olarizatio n pro- 
files. It has been shown by S anchez Almeida et all (| 19961 
§ 4.2) that, under the MISMA hypothesis, three is the min- 
imum number of components with constant magnetic field 
required to reproduce the asymmetries of the Stokes profiles 
observed in the network, a nd three seem to suffice to repro - 
duce all quiet Sun profiles ([Sanchez Almeida fc Litesll2000f ). 
Such components should be considered as a schematic rep- 
resentation of the average properties of the atmosphere (see 
Sanchez Almeida fc Lites 2000, § 5). Each component of 
the model is characterized by the variation with height of 
the thermodynamical properties, plasma motions along the 
magnetic field lines, occupation fraction, and magnetic field 
strengths. These variations are forced to conserve magnetic 
flux and mass flows. The temperature stratification is im- 
posed to be the same in the three optically thin compo- 
nents, since we expect an efficient radiative thermal ex- 



4 



Viticchie et al.: MISMA interpretation of HINODE SOT/SP data 



change among them. Moreover, the lateral pressure bal- 
ance links the properties of the components at every height; 
they must have the same total pressure defined as the sum 
of the gas pressure plus the magnetic pressure. The me- 
chanical balance couples the magnetic field and the ther- 
modynamics in the model atmosphere. Magnetic field in- 
clination and azimuth are constant with height. Finally, 
an unpolarized stray-light contribution is considered (for a 
complet e description of the inversion code and hypothe- 
ses, s ee ISanchez Almeidal Il997l : ISanchez Almeida fc LitesI 
I2QQQD . The inversion of the some 400 wavelengths defining 
each set of Stokes profiles employs only 20 free parameters, 
which still are twice the number of fre e parameters of ME 
inversions (|Orozco Suarez et al.ll20073 ). 

The adopted inversion strategy is the following. Stokes 
/ and V are inverted in those pixels presenting \V\ ampli- 
tude above 4.5 x cjy at least for one wavelength in one 
of the two Fel lines. If \Q\ or \U\ amplitude is above 
4.5 X cfq^Ui at least for one wavelength in one of the 
two Fel lines, a full Stokes inversion is also performed. 
According to these criteria about 29 % of the subfield has 
been inverted, with 2.3 % corresponding to full Stokes in- 
versions. The 4.5 X (JQ^uy thresholds have been chosen 
so as to be the same as that of the ME inversion analy- 
ses already performed on our dataset (O rozco Suarez et al.l 
Hob 7a; Asensio Ramos 2009). Each profile has been inverted 
using as initial ^ues s all the elev e n clas ses of MISMA mod- 
els of .Sanchez Almeida fc LitesI (|2000l ) . Among the eleven 
inversions, we select the one having the smallest deviation 
between the observed and the computed profiles. Note that 
when only Stokes / and V are inverted, we assume the mag- 
netic field to be longitudinal for both the magnetic compo- 
nents (i.e., the inclination is either zero or 180°). As we 
argue in Appendix [A| this assumption influences the frac- 
tion of stray light inferred from the inversion, but not the 
magnetic field strength. 

3.1. Inversion examples 

In this section we present examples of inversions. They have 
been selected among all the 11600 inverted profiles, and 
they are representative of the goodness of the fit. 

Fig. [3] shows three examples of pixels where the inver- 
sion code retrieves two opposite polarity fields coexisting in 
the resolution element. They have been picked in different 
regions of the subfield, i.e., around the negative network 
patch on the top-left corner of Fig. [2l and in IN regions. 
The profiles show how the two polarities can be either eas- 
ily detectable, when quasi symmetric Stokes V profiles are 
measured, or almost hidden, when one of the two lobes of 
Stokes V almost vanishes. Two of the three selected cases 
coincide with a spatially-resolved large-scale change of po- 
larity, i.e., the inverted pixels correspond to frontiers be- 
tween positive and negative polarities in the COG mag- 
netogram. In addition, the third example in Fig. [3] shows 
that mixed polarities can be found associated with polar- 
ized pixels that are not in frontiers between opposite polar- 
ities. Note that the three cases reported in Fig. [3] show no 
evidence of linear polarization signals above the 4.5 x cfq^u 
threshold. Even below the imposed threshold, the linear po- 
larization profiles do not present any shape that could be 
interpreted through an inversion procedure with reliability. 

Figure [4] contains three examples of profiles observed far 
from network regions. They have been chosen to be char- 
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Fig. 3. Examples of MISMA inversion of HINODE 
SOT/SP Stokes profiles corresponding to mixed polarities 
in the resolution element. We show profiles in three different 
pixels representative of a weakly polarized IN region (upper 
panels), and frontier pixels between regions with opposite 
polarity in the IN (central panels) and network (lower pan- 
els). The results are organized in two sets of panels; each 
of these sets contains an upper row with images that put 
the pixel into context, and a lower row with the four Stokes 
profiles (as labeled). The upper row shows I x I {I c:^ 6") 
maps of the COG magnetogram saturated at ±50 G (left), 
the continuum intensity (central), and the net-circular po- 
larization saturated at 0.3 (right). The red square indicates 
the position from which the profiles are taken. Lower rows: 
observed Stokes profiles (symbols) and best-fitting Stokes 
profiles (solid lines). The vertical dotted lines mark the lab- 
oratory wavelengths of the two lines. 



acteristic of the typical profiles found in the IN as observed 
by HINODE SOT/SP. All Stokes V profiles present asym- 
metries, for which the inversion works in a really satisfac- 
tory way. The central panels of Figs. [3] and H] correspond 
to adjacent pixels next to an apparent neutral line in COG 
magnetogram (i.e., where the magnetic field changes sign; 
compare their position on the magnetogram). They have 
been selected to show how the MISMA inversion discerns 
between a single polarity model (Fig. |4]) and a mixed polar- 
ity model (Fig. [3]). The Stokes V profile in the central panel 
of Fig. [3] presents a dominant negative blue lobe which in- 
dicates the presence of negative fields. Beside this, a small 
negative red lobe is also detected by the inversion. The 
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Fig. 4. Examples of MISMA inversion of HINODE 
SOT/SP Stokes profiles corresponding to a single polar- 
ity in the resolution element. For details on the layout of 
the figure, see the caption of Fig. [3l 





Fig. 5. Examples of full Stokes MISMA inversion of 
HINODE SOT/SP Stokes profiles. For details on the layout 
of the figure, see the caption of Fig. [3l 



latter denotes the presence of positive fields. The MISMA 
code correctly interprets such a profile as emerging from 
a pixel in which opposite polarities coexist. Such an inter- 
pretation is consistent with the COG map, showing how 
the pixel from which the profile is taken lies in the frontier 
between positive and negative field patches. The central 
panel of Fig. S] shows the Stokes V profile of a pixel next 
to the previous one in the direction of the positive polarity 
concentration. It is very different from its neighbor - the 
profile in Fig. H] is dominated by a negative red lobe and 
presents a small positive blue lobe. The MISMA code in- 
terprets such a profile as emerging from a pixel in which 
only positive fields are present. In this case the inversion is 
also consistent with the COG magnetogram, which shows 
the pixel lying on a positive polarity patch. 

Figure [5] presents two examples of full Stokes inversion, 
i.e., inversions including Q and U. Such examples illustrate 
the linear polarization signals we considered to be invert- 
ible. The upper panel of Fig.[5]shows the inversion of a pixel 
very close to the example in the upper panel of Fig. [3l In 
this case the linear polarization signal is strong enough to 
be analyzed and the MISMA code succeeds in the inversion. 
It is important to notice that in such a pixel, similarly to 
what found in the upper panel of Fig. [3l polarization pro- 
files are still interpreted by the code as emerging from a 
mixed polarity pixel. In the lower panel of Fig. [5] a full in- 
version of an IN pixel is represented. In this case, even if 
the pixel is almost on the frontier between opposite polarity 
regions, a single polarity is measured. In the two cases the 
magnetic field strengths at the base of the photosphere are 
in the kC regime. 

The examples discussed here illustrate not only the 
goodness of the fits but also the soundness of the MISMA 
interpretation of HINODE SOT/SP measurements. Such 
measurements are often characterized by important asym- 
metries in Stokes V profiles; Figs. [3] and |4] show three ex- 
amples of Stokes V profiles whose NCP > 0.3. From the 
maps in Fig. [2j we notice that such values for the NCP 
are very common in the selected 29.52'' x 31.70'' subfield 
and, by extrapolation, they should be very common in the 
full FOV as well. The common presence of large asymme- 
tries demands a refined inversion method to interpret quiet 
Sun SOT/SP profiles, such as the MISMA inversion we are 
employing. Detail on the percentage of asymmetric profiles 
are reported in § SI 



4. Results 

The inversion code succeeded in the inversion of 11600 pro- 
files which represent 29% of the selected subfield. The total 
time needed to perform such an analysis is about two days 
when the inversion analysis is split into eight IDL jobs run- 
ning in parallel 0. Among the 11600 inverted profiles, 925 
(::^ 8 % of the total number) show invertible linear po- 
larization signals. These profiles have been found mainly 
in correspondence with patches of strong polarization (see 
the upper panel in Fig. [5]). Here we present the results ob- 
tained from the inversion of Stokes / and V profiles alone. 
Full Stokes inversions using modified MISMA models hav- 



^ The inversion of each pixel is independent so the analysis 
of the invertible domain can be distributed over different CPUs; 
the jobs run on two Xeon E5410 Quad-core processors with 8 GB 
of shared ram. 
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Fig. 6. Results of the inversion at the constant reference height, corresponding to the base of the photosphere. The differ- 
ent panels contain: classes of initial model MISMAs used in the inversion (upper- left panel), field strength for the minor 
component (upper-central panel), occupation fraction for the minor component (upper-right panel), COG magnetogram 
saturated at ±50 G (lower- left panel), field strength for the major component (lower-central panel), occupation fraction 
of the major component (lower-right panel). The contours on the COG magnetogram represent the network regions as 
identified by our automatic procedure (red), and the pixels with mixed polarities in the resolution element (blue). Black 
pixels in the MISMA class, field strength, and occupation fraction maps correspond to regions that are not inverted (see 
Fig. [2]). The bars on top represent the color palettes adopted for the MISMA initial models (on the left), and for the two 
pairs of images of field strength (on the center) and occupation fraction (on the right) 



ing inclination and azimuth as free parameters have been 
performed with success (e.g., Fig.[5|), but they are less rep- 
resentative from a statistical point of view, and they pro- 
vide similar results. In fact, magnetic field strengths from 
few hC to kC fields are measured in fully inverted pixels, 
with a slightly higher probability for kC fields. We focus 
our analysis on the statistical properties of the magnetic 
field strength and occupation fraction at the photosphere 
as retrieved by the inversion code (the occupation fraction 
is just the volume filling factor). Assuming the model atmo- 
spheres to be in lateral pressure balance, each value of total 
pressure defines the same geometrical height in all model 
MISMAs. We decided to take as reference height in the 
atmosphere the one corresponding to the base of the pho- 
tosphere in ty pical ID quiet Sun model atmospheres (e.g., 
iMaltbv et al.| [l986). The photospheric lower boundary in 
such models is defined as the height where the continuum 
optical depth at 500 nm equals one, which corresponds to 
a pressure of about 1.3 x 10^ dyn cm~^. Unless otherwise 
is mentioned explicitly, all the parameters discussed here- 
after refer to this reference height or, equivalent ly, to this 
reference total pressure. 

Fig. [6] shows six maps summarizing the inversion at 
the reference height, together with the COG magnetogram, 
used here to show the context. The red contours on the 
COG magnetogram separate network and IN regions in the 
selected subfield. Network patches have been identified us- 
ing an algorithm which takes into account two properties 
of the IN and network regions. First, the IN covers most 
of the solar photosp here at any time (more than 90 % for 
iHarvev- Angle! [19931 which is the figure we use in the cal- 



culations). Second, network patches are magnetic concen- 
trations showing a spatial coherence. The network must 
be identified not just as a strongly polarized region; it is 
expected to be also spatially extended. The procedure to 
define the network works as follows: (1) it automatically 
finds the threshold on the total polarization that makes the 
patches of large signal to cover a few percent of the FOV 
These signals are tentatively identified as cores of 
network elements. (2) It checks the spatial extension of the 
network patch candidates so as to exclude connected struc- 
tures composed by less than ten pixels. Once the network 
patches have been identified, the procedure dilates them 
using a 10 X 10 pixel square kernel. The network patches 
thus selected cover 11 % of the subfield. 

At the location of the network patches the map show- 
ing the initial guess models for the inversion of each pixel 
reveals a significant degree of local coherence (Fig. [6l top 
left panel). On the contrary, the initial model MISMA fluc- 
tuates in the IN on scales smaller than 1" . Considering net- 
work and IN altogether, the percentages of different types 
of best initial guess models are class ~ 5%, class 1 18%, 
class 2 c:^ 16%, class 3 :^ 10%, class 4 c:^ 11%, class 5 c:^ 0%, 
class 6 ^ 7%, class 7 12%, class 8 14%, class 9 5%, 
and class 10 0% (see Sanchez Almeida & Lites 2000 
for a detailed description of the classes). The initial model 
MISMAs can be used as a proxy for the kind of asymme- 
tries present in HINODE spectra. We find that 35 % of the 
inverted profiles belong to classes 4, 6, 7, 9. These are the 
classes where Stokes V presents large asymmetries, similar 
to those shown in Figs. [3] and [H 
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The central panels of Fig. [6] show the magnetic field 
strengths for the two magnetized components of the model 
MISMAs. The results have been organized so that the ma- 
jor and the minor components are considered separatel}|3. 
On the one hand, the major component presents occupation 
fractions (equivalent to volume filling factors) going from 
~ 10~^ up to > 10"-^, and magnetic field strengths between 
G and 1.8 kG. On the other hand, the minor compo- 
nent presents filling factor almost always ^ 10~^ and field 
strengths ^ 1 kG. Only in the network patches the minor 
component presents values for the filling factor comparable 
with the values of the major component. 

The statistical properties of the field strength for both 
the major and minor components are reported in Fig. [71 
It gives the fraction of analyzed photosphere covered by a 
given field strength Pz{B)^ i.e., except for a scaling factor, 
it gives the sum of all occupation fractions for each mag- 
netic field strength B. The histograms have been defined 
by sampling the interval — 2 kG with 20 bins of 100 G. 
They show the probability of finding a given field strength 
in our quiet Sun inversion. The distribution corresponding 
to the whole subfield presents a maximum at o:^ 1.6 kG, 
and then an extended tail for smaller field strengths. The 
statistics for the IN (Fig. [71 top, thick solid line) also shows 
kG fields, with a secondary maximum at the equipartition 
value of 5 hG. Note the reduction of the kG peak with 
respect to the full subfield distribution (Fig. [71 top, thin 
solid line). It is clear that the hG fields are almost com- 
pletely in the IN, while the network is dominated by kG 
fields with little contribution for fields in the hG regime. 
On the other hand, the statistics of the minor component 
is completely dominated by strong kG fields with a maxi- 
mum just at 1.8 kG and almost null hG contribution; this 
is valid for the whole subfield, IN, and network. 

The magnetic field strength distribution in Fig. [8] is 
calculated for a height in the atmosphere of approxima- 
tively 150 km, i.e., a height representative of the formation 
region of the core of the Fel lines. Similarly to what is 
done for the definition of the reference height of the pho- 
tosphere, the heights are selected as those having a total 
pressure of 4 x 10^ dyn cm~^ (approximatively at 150 km 
in the pressure stra tification of the mean photosphere; e.g. 
'Maltby et al. 1986). The main difference between the two 
statistics in Figs. [71 and [8] is the large reduction of kG fields 
at 150 km. The magnetic structures are treated as thin fiux- 
tubes by the inversion code, therefore, the magnetic field 
lines must fan out with height to maintain the mechanical 
equilibrium, which implies a drop in field streng th, and a n 
expansion of the magnetic structures (e.g., SpruitI [TosTI ) . 
The two factors conspire to produce the observed increase 
of hG fields with height (cf. Fig.[71and Fig.[8|). For example, 
the volume occupied by magnetic fields smaller than 1 kG 
at 150 km is some 8 times larger than their volume at the 
base of the photosphere. 

A very interesting result, unique to MISMA inversions, 
is the detection of a large number of mixed polarity pixels 
(three examples are reported in Fig. [3]). These are identified 
by the blue contours on the COG magnetogram of Fig. [6l 
Pixels presenting mixed polarities are very common in IN, 
namely about 25 % of the total number of inverted profiles. 



^ In each pixel, the major and the minor components are 
the magnetized components having the largest and the small- 
est mass, respectively. 
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On the contrary, network patches tend to be interpreted by 
the code as unipolar regions. We note that mixed polari- 
ties show up almost always either in the frontiers between 
strong signals of opposite polarity, or in low polarization re- 
gions, close to the 4.5 x ay level. We also note the presence 
of mixed polarities surrounding the large negative polarity 
network patch in the selected subfield (see the upper left 
corner in the COG magnetogram in Fig. [6]). They are al- 
ways found where positive small concentrations are close to 
the network patch. Fig. [71 also includes the statistical prop- 
erties of the mixed polarity pixels separately. Contrarily 
to the other statistics described above, the distribution for 
mixed polarity pixels is not dominated by kG fields: hG 
and kG fields are found almost with the same occupation 
fraction in both IN and network regions, and for both mi- 
nor and major components. The distribution presents an 
absolute maximum at about 1.6 kG and a secondary one 
at about 500 G. 

The fraction of analyzed photosphere occupied by mag- 
netic fields a can be computed as the integral of the distri- 
butions in Fig. [71 i.e., 

a= / Pz{B)dB. (1) 
Jo 

Considering the full subfield, the value of a for the ma- 
jor and the minor components is 3.1 % and 1.4 %, re- 
spectively. In other words, 95.5 % of the analyzed photo- 
sphere is field free or, more precisely, it has field strengths 
much smaller than the ones inferred from inversion, so 
that the inversion code cannot distinguish them from zero. 
Other fractions corresponding to the various components 
in the FOV are listed in Table [H The first moment 
of the magnetic field strength distributions also provides 
an estimate of the unsigned fiu x density (B) (see, e.g., 
iDommguez Cerdeha et al.ll2006b[ ), 

POO 

{B)= / BPz{B)dB, (2) 
and of the average field strength, 

B={B)/a. (3) 

The values of {B) and B for the various components in the 
FOV are listed in Table [H As the table shows, {B) = 66 G 
considering network and IN all together. This fiux den- 
sity is so large partly because it refiects the fraction of 
the 29. 52'^ x 31.70^' subfield having the largest polariza- 
tion signals. Such fraction / is about 0.29, therefore, if one 
considers the full subfield, the fraction of photosphere oc- 
cupied by the magnetic fields inferred from our inversion is 
a/, and the unsigned fiux density {B)f. The values of the 
filling factor and unsigned flux densities corresponding to 
this other normalization are also included in Table [1] within 
parentheses. Then the unsigned flux density decreases to 
about 19 G. This still large value is put into context in §[5l 
but it is important to realize that it implicitly assumes the 
(1 — /) non-inverted subfleld to have no magnetic flelds. 

We have also tried a crude separation between gran- 
ules and inter granules. A simple thresholding criterion was 
used, so that pixels brighter than the mean intensity are 
granules, and vice- versa. We flnd that 75 % of the magne- 
tized plasma is in intergranules, and this fraction increases 
with increasing fleld strength - 85 % of the plasma having 
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B > l.b kG is in intergranules. These figures are meant to 
be rough estimates since our simple criterion is insufficient 
for an accurate separation between granules and intergran- 
ules. 
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Fig. 7. Fraction of analyzed photosphere covered by a 
given strength. Field strength distributions for the major 
component (upper panel): IN (black thick line), network 
(red line), mixed polarity regions (blue line), and the full 
29.52'' X 31.70'' subfield (black thin line). Field strength 
distributions for the minor component (lower panel). 
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Fig. 8. Fraction of analyzed photosphere covered by a given 
strength at 150 km height in MISMA models. For details 
on the layout of the figure, see the caption of Fig. [71 

5. Discussion 

The ubiquity of Stokes profiles with asymmetries implies 
that even HINODE SOT/SP resolution is insufficient to 



resolve the structure of the q uiet Sun magnetic fields (see 
ISanchez A lmeida et al.l Il996f ). This unresolved structure 
has been neglected in the analyses carried out so far, and 
we present the first attempt to account for it when inter- 
preting HINODE spectra. We em ploy the MISMA inversion 
code bv lSanchez Almeidal (|1997[ ). which produces very good 
fits. The MISMA description is able to reproduce the whole 
range of profile shapes observed by HINODE SOT/SP in 
the quiet Sun, which are often extremely asymmetric. This 
was known for Advanced Stokes Polarimeter (ASP) spectra 
(Sa nchez Almeida fc LitesI l2000f ). but HINODE SOT/SP 
has three times the spatial resolution of ASP (i.e., 0.3" 
against 1"). 

The percentages of MISMA classes used as ini- 
tial guess are different from the ones reported in 
ISanchez Almeida fc LitesI (j2000[ ). This can be due either 
to the inversion strategy (e.g., the imposed threshold to se- 
lect invertible data) or to a general modification of Stokes 
profile shapes due to the higher spatial resolution. 

As pointed out in ISanchez Almeida fc Lite's! (|20 00V the 
major component is the dominant source of polarization of 
the model since it contains the largest mass. For this rea- 
son, the statistics of the major component will be consid- 
ered as the reference statistics derived from our analysis. 
The field strength distribution for the major component 
sh ows an important difference wit h respect to the results 
in ISanchez Almeida fc LitesI (|2000f ): in fact, our analysis is 
able to reveal much more hG fields than the analysis of 
the ASP data. These are almost completely in the IN, and 
a large part of them are found in mixed polarity pixels. 
Moreover, the statistics on the latter shows that hG and 
kG fields are measured with the same probability. The IN 
kG magnetic fields are more i mportant in our invers i on than 
the inversions carr ied out bv lOrozco Suarez et aP (l2007a[ ) 
and ^ Asensio RamosI (|2009 ) . Part of this difference is cer- 
tainly due to their use of ME atmospheres, unable to repro- 
duce the important line asymmetries characteristic of the 
quiet Sun Stokes profiles. Such difficulty may be secondary 
when dealing with slightly asymmetrical profiles, but the 
field strengths are inferred from the Stokes V shapes, and 
any field strength extracted by fitting anti-symmetric ME 
profiles to profiles like those in Figs. [3] and |4] is open to 
question. Another part of the difference can come from 
the large drop of magnetic fiel d strength with height in 
the atmosphere. As discussed in lDommguez Cerdeha et al.l 
(f200 6b) and in § SI a reduction of the magnetic pressure is 
required to maintain the mechanical balance between the 
magnetic structures and the field free atmosphere when the 
gas pressure decreases exponentially with height. Such a de- 
crease in field strength is intrinsic to the MISMA models, 
which force their component s to be in horizontal a nd ver - 
tical mechanical equilibrium. l Orozco Suarez et al.l (|2007b[ ) 
showed that the ME inversion of Fel 630 nm gives us in- 
formation on the solar atmosphere at the height of forma- 
tion of the lines, which is higher than the reference height 
adopted in this work. When the magnetic field strength 
statistics is calculated at a height of 150 km, then kG fields 
are almost absent (Fig. [8]), which partly reconciles our re- 
sults with those inferred from ME inversions. However, the 
shap e of our distribution still d iffers from the one reported 
in Orozco Suarez et al.l (|2007al ) in which the most probable 
field strength is :^ 100 G. 

Three additional comments regarding the field strength 
distribution are in order. First, kG fields are needed to 
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Table 1. Summary of magnetic properties at the base of the photosphere derived from the MISMA inversion 



Portion of photosphere 


Field Streni 


5th B [G] 




Flux Density {B) [G] 






Filling Factor 


a 






major 


minor 


major minor 


both 


major 


minor 




both 


internetwork (IN) 


1258 


1521 


20 


(5.8) 12 (3.5) 32 


(9.3) 


1.6 (0.5) 


0.7 (0.2) 


2.3 


(0.7) 


mixed polarity regions 


1036 


1041 


5 


(1.6) 1 (0.3) 6 


(1.9) 


0.5 (0.1) 


0.1 (0.0) 


0.6 


(0.1) 


network 


1513 


1644 


23 


(6.7) 11 (3.2) 34 


(9.9) 


1.5 (0.4) 


0.7 (0.2) 


2.2 


(0.6) 


full subfield 


1380 


1578 


43 


(12.5) 23 (6.7) 66 


(19.2) 


3.1 (0.9) 


1.4 (0.4) 


4.5 


(1.3) 



Note 1. The four rows are defined from the four distributions in Fig.d which are normalized to the analyzed portion of photosphere, 
i.e., 29% of the 29.52^^ x 31.70^^ subfield. The values normalized to the whole subfield area are reported within parentheses. 



explain the presence of G-band bright points in the quiet 
Sun intergranular lanes (e.g., Sanchez Almeida et al. 
2004 Ide Wiin et all l2005t iBovelet fc Wiehri 120081 : 
Sanchez Almeida et al.l |2010[ ) and, therefore, it is a 
sign of consistency that they show up in spectro- 
polarimetric observations. The occupation fraction of kG 
fields inferred from the inversion is of the order of a few 
percent, which is in good agreement with the 0.5 — 2 % 
G-band bright point su r face c o verage recently ni e asured 



(ISanchez Almeida et all I2QQ4I : iBovelet fc Wiehrl I2QQ8I : 



ISanchez Alm eida et al." "2Q1QV Moreover, the large polar- 
ization signals that we select for inversion appear preferen- 
tially in the intergranular lanes (|Dommguez Cerdefia et al.l 
l2QQ3al : iLites et all l2QQ8f ) where the BPs are known to 
reside. Second, our results do not exclude the presence of 
weak fields (i.e., fields producing weak polarization signals 
in the Fel visible lines). Our inversion retrieves the field 
strength of only 4.5 % of the photospheric plasma (see §[4|). 
The field strength of the rest remains unconstrained. This 
part is almost certainly magnetic, but with properties that 
remain elusive to our observation because of its complex 
magnetic topology that cancels Z eeman signals, because 
of the low fi eld strengths, or both. Sanchez Almeida et al.l 
([2003) and Dommgue z Cerdeha et a l. (2006a') show how 
part of these hidden fields can be revealed through MISMA 
inversion by analyzing Fel infra-red lines. Even weaker 
fields can be inferred and diagnosed via Hanle effect (e.g., 
iFaurobert-SchollI Il994l : iBommier etld] I2QQ5). In order to 
avoid overlooking the weak field component, the quiet 
Sun regions need to be analyzed using measurements in 
both visible and infra-red spectral ranges, and produced 
by both Zeeman and Hanle effects. How to consistently 
combine these measurements remains an open issue, but 
such combination is certainly the way to go (se e the 
first attempt by iDommguez Cerdefia et all l2QQ6b[ ). The 
third comment on the field strength distribution refers to 
discarding biases that may art ificially produce kG fields. 
iBellot Rubio fc ColladosI ((2QQ3D show how large noise m 
the Stokes V profiles of the Fel 630 nm lines may fake 
kG fields. However, this effect is insignificant above a 
signal-to-noise ratio threshold that our typical spectra 
exceed. Moreover, in addition to kG, noise also produces 
false weak dG, which we do not see. Another potential 
source of false kG fields could be (polarized) stray-light 
from the network to the IN. Although we cannot discard 
some influence close to the network patches, stray-light 
cannot account for the bulk of the kG fields in the IN. 
Stray-light is at most 10% as me asured for the broa dband 
filter imager aboard HINODE (|Wedemever-BQhml l2QQ8t) . 
It could account for 10% of the network magnetic fiux 
artificially appearing in the IN, but this represents only a 



few G, thus unable to explain the some 30 G inferred for 
the IN (see Table [I]). 

The mixed polarity regions revealed by the MISMA 
analysis are about 25 % of the total number of in- 
verted profiles, similar t o the percentage reported in 
Sanchez Almeida fc Lite s (I2QQQD f or ASP data (see also 

I So cas- Navarro fc Sanchez Almeida! l2QQ2f ). Note that the 
percentage of mixed polarities has not changed despite the 
increase of angular resolution (from 1^' to 0.3'^), and despite 
the fact that the polarization signals analyzed here are gen- 
erally larger (we use a noise threshold 1.5 times larger than 
the one used on ASP spectra). The really significant differ- 
ence between the two datasets is t he fraction of solar surface 
covered by polarization signals; ISanchez Almeida fc LitesI 
(2000) invert only 15 % of the surface, whereas HINODE 
spectra allow us to almost double this fraction. The in- 
crease of angular resolution has certainly resolved some of 
the mixed polarities detected at ASP angular resolution, 
however, at the same time, the newly revealed Stokes V 
profiles uncover mixed polarities at even finer spatial scales. 
This observation is consistent with the numerical models of 
solar mag neto- convection and/or tu rbulent dynamos (e.g., 
ICattaned[T999l : IVoder fc Schusslerll2QQ7[ ). They predict an 
intricate pattern of highly intermittent fields varying over 
very smah spatial scales down to the diffusion length-scales. 

The average flux density we infer from the inverted 
profiles is of the order of 19 G when the area of the 
full subfield is considered for normalization (see Tabled]). 
Thi s value is almost twice a s larg e as the one inferred 
by ISanchez Almeida fc Lit esI (l20Q0|, 1" resolution), and 
it is also significantly larger than the ones obtained 
from ME inversions of qui et Sun HINODE spectra by 
lOrozco Suarez et al.l (l2007al - 9.5 G f or the fuh FOV) 
and using the magnetograph equation by Lites et al. (2008', 

II G). However, our estimate is comparable with other 
HINODE based estimates of the unsigned flux of the quiet 
Sun (e.g., Jin et al. 2009 get 28 G). Such a difference could 
be explained taking into account that ME inversions are 
blind to sub-pixel magnetic structuring. Note that the av- 
erage flux density from the major component (i.e., the dom- 
inant source of polarization ) is i n good agr e emen t with 
lOrozco Suarez et al.l (|2QQ7aD and iLites et al.l (|2QQ8l ) - we 
get 12.5 G, see Table [H The small excess with respect to 
their results is compatible with the 1.6 G flux contribu- 
tion from mixed polarities (see Table [T]). The bulk of the 
difference seems to reside in the 6.7 G provided by the mi- 
nor component. Its contribution to the flux density is more 
important than the modification it produces on the polar- 
ization signals (even though it is needed to reproduce the 
line asymmetries). 



10 



Viticchie et al.: MISMA interpretation of HINODE SOT/SP data 



6. Conclusions 

We present the first inversion of HINODE SOT/SP Stokes 
profiles tfiat accounts for tfie asymmetries of tfie profiles 
in the quiet Sun. The analysis is carried out under the 
MISMA hypotheses which allows to reproduce the dif- 
ferent profile shapes observed in the quiet Sun. We fol- 
low the approach already successfully exploited to de- 
scribe the asymmetries observed with 1" angular resolution 
(jSdnchez Almeida & Lites 2000). 
The main results are as follows: 

1. We inverted 11600 sets of Stokes / and V profiles rep- 
resentative of the quiet Sun internetwork (IN) and net- 
work at disk center. The inversion code is able to repro- 
duce, in a satisfactory way, the whole variety of asym- 
metries revealed by HINODE SOT/SP. 

2. The MISMA code is also able to reproduce linear po- 
larization measurements when full Stokes inversions are 
performed. Here we report a few examples. 

3. The existence of asymmetries is certainly not negligi- 
ble. Some 35 % of the analyzed profiles present large 
as ymmetries, according to the roug h classification used 
in Sanchez Almeida fc Li tes' (2000). 

4. 25% of the analyzed profiles present asymmetries that 
are interpreted by the MISMA code as due to mixed 
polarities in the resolution element. These pixels are 
found to be located either in transition regions between 
patches of opposite polarity, or in pixels presenting weak 
polarization signals. 

5. The magnetic plasma whose properties the inversion 
code constrains represents only 4.5 % of the photo- 
spheric plasma. The rest remains unconstrained by our 
analysis. 

6. The statistics of the detected magnetic field strength at 
the formation height of the continuum is dominated by 
strong kG fields, both in the network and the IN. The 
later, however, presents an extended tail over the whole 
hG domain. 

7. At the height of 150 km (representative of the forma- 
tion heights of the Fe I visible lines) very little kG fields 
are measured. This result narrows down the gap be- 
tween our field strength distribution and those inferred 
from ME invers i ons bylOrozco Suarez et al.l (|2007af ) and 
lAsensio RamosI (|2009[ ), and it follows directly from the 
decrease of the field strength in response to the decrease 
of the gas pressure with height in the photosphere. 

8. The average flux density derived from the inverted pix- 
els is 66 G. If one considers the full analyzed subfield, 
it becomes 19 G. The same figures for the IN become 
32 G and 9.3 G, respectively. These values are signifi- 
cantly larger than the ones obtained from ME inversions 
of quie t Sun HINODE spectra by Orozco Suarez et al. 
(|2007al . - 9.5 G for t he full FOV ) , and using; the magne- 
tograph equation bv llLites et al.l (|2008L 11 G). However, 
our estimate is comparable with other HINODE based 
estimates of the unsigned flux of the quiet Sun (e.g., Jin 
et al. 2009 get 28 G). 

As a general concluding remark, the analysis here pre- 
sented reveals the importance of a complete interpretation 
of the shape of Stokes V profiles through refined inver- 
sion codes. This work should be considered as the first step 
in the interpretation of HINODE SOT/SP profile shapes. 
Because of the wealth of information contained in the asym- 



metries, such a topic deserves attention from the solar com- 
munity. Different inversion codes and hypotheses can be 
adopted to interpret such asymmetries, and they will prob- 
ably be able to reproduce the line shapes as well (e.g., a sys- 
tema tic variation of magnetic field and v elocity along the 
LOS: lRuiz Cobo fc del Toro Iniestal [19921 ). The uniqueness 
of the interpretation cannot be assessed unless different al- 
ternatives agree. 
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Fig. A.l. Two sets of synthetic Stokes / and V profiles 
of Fel 630.25 nm obtained with different magnetic field 
inclinations but the same product (1 — fsi) cos^, with fsi 
the fraction of stray-light. The solid and the dotted lines 
correspond to = 0°and 60°, respectively. The two sets 
cannot be distinguished within typical HINODE noise. 
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Appendix A: Lack of crosstalk between magnetic 
field strength and magnetic field inclination 

When Q and U are too small to be used, the inversions assume 
the magnetic field to be longitudinal. It is important to realize that 
such assumption does not influence our magnetic fleld strength di- 
agnostics. When the polarizations signals are weak (our case), and 
the inclination is constant (our assumption), then Stokes V scales 
with the cosine of the inclination in dependently of the fleld strength 
(jSanchez Almeida Truiillo BueiiQi ri999, § 3.1). The information on 
the fleld strength is coded in the shape of Stokes whereas the mag- 
netic fleld inclination modifles the scaling. Such scaling is indepen- 
dently obtained by the inversion code disguised as stray-light, there- 
fore the inferred magnetic fleld strength and inclination are uncou- 
pled. All quiet Sun inversions assume the pixels to be partly covered 
by magnetic flelds, and the corresponding fllling factor is determined 
as a free parameter of the inversion (stray-light factor). An error on 
the magnetic fleld inclination affects the stray-light fact or, b ut it does 
not influence the magnetic fleld determination. Figure [A. II illustrate 
the argument. It shows Stokes / and V synthesized in two atmo- 
spheres that differ in stray-light factor fgi and cosine of inclination 
cos^, but they have the same product cos^(l — fsi). The two pairs 
of Stokes proflles are indistinguishable within the noise of HINODE 
observations. 



